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Description 

FIELD OF THE INVENTION 

5 [0001 ] The present invention is referred to the field of radiomobile systems and more precisely to a variable complexity 
space-time channel estimation and receiver for mobile communication systems. 

BACKGROUND ART 

10 [0002] Main problem of digital receivers is that to correctly detect the transmitted sequences of data, that because 
some electromagnetic phenomena having casual nature distort the waveforms of the transmitted signal making them 
not more immediately recognizable to the receiver. Noise always present on the channel and interferences worsen the 
received signal furthermore. So the receiver shall eliminate from the received signal the effect of noise and interferences 
to achieve a correct data detection. In the field of the invention the performance of these tasks needs the preliminary 

15 estimation of the channel response in concomitance with the arrival of a training sequence, known a priori by the 
receiver. As far as concerns thermal and floor noise, the same is considered spectrally white, gaussian in the amplitudes 
(AWGN) and temporally uncorrelated to the useful signal. Data detection can avail, for example, of a linear receiver 
which attempts to de-correlate the transmitted data from the equalized reception, or make provision for a Viterbi proc- 
essor which performs a Maximum Likelihood Sequence Estimate (MLSE). In such a case it doesn't attempt to equalize 

20 the channel in strict sense, but rather uses the knowledge of the channel pulse response (get from the training sequence 
estimation) to find the data sequence transmitted with the maximum probability. To summarize in a practical case 
concerning the well known GSM, the receiver firstly estimate the channel response during the reception of the midam- 
ble, then the estimate is frozen for all the burst duration and used for the detection of the remaining semi-bursts of data 
at the two sides of the midamble. This way of operation is strictly valid for stationary channels and low mobile speed. 

25 A more realistic approach takes into account the variability of the channel and attempts to extend the channel estimation 
even outside the midamble in order to pursue said variability and compensate the mobile speed. 
[0003] As known, a main cause of distortion is the multipath fading due to the scattering of the electromagnetic waves 
from atmospheric diff users statistically spread along the radio propagation. Flat or selective fading are generated con- 
sequently. The first involves the whole band and is due to small delays along the various paths, being the delays each 

30 other correlated and consequently assumed like a phase offset. Selective fading is due to uncorrelated delays com- 
parable to the inverse of the band width; the observed depths inside duration intervals obey to a Rayleigh distribution. 
The most popular countermeasures against the Rayleigh fading are diagonal interleaving of the code blocks on more 
transmission bursts, active against selective fading, and the space diversity reception by using at least two antennas 
spaced more than X/2 t assuming uncorrelated fading on the two paths. 

35 [0004] Recently the opportunity to exploit directive antenna arrays in the Base Station (BTS) of radiomobile systems 
has been taken into consideration and a certain number of patent applications have been filed on the argument. The 
directive array theory, in particular the Uniform Linear Arrays (ULA) was originally developed in radar and sonar military 
applications to the precise aim of localizing the positions of possible targets. A potential target reflects along some 
directions either a radio or acoustic beam appositely transmitted and its localization involves the estimate of the Di- 
rection Of Arrival (DOA) of the reflected waves. The DO A estimation has been also proposed in the radiomobile ambit 
for associating a temporal channel to each estimated direction of arrival of a useful signal at the BTS side. Such a 
channel, supposed be symmetric along the two directions, allows a double advantage: saving transmission power both 
from the Mobile Stations and the BTS because the power is prevalently focussed towards the receivers, and the re- 
jection of cochannel interferents from the BTS only by means of an opportune reception beamforming. Other inventions 

45 also based on the DOA estimate (i.e. in the name of ARRAYCOM) allow to implement a Space Division Multiple Access 
(SDMA) technique, by which more than one mobile station MS are allocated on the same carrier, and the MSs are 
discriminate each other from the respective DOA. Nevertheless the DOA estimate is very complicate and only profitably 
in case of exploitation of the SDMA technique, while it appears cumbersome elsewhere, i.e. at the only aim of estimating 
the channel response. 

so [0005] For the last aim the German patent application DE-A-1 963941 4, inventor: RADEMACHER L. discloses a BTS 
receiver able to jointly estimate the space-time channel without recurring to a separate DOA estimate. The claim 1 of 
this citation recites textually: 

"Method for the parametrization of a receiving station (BS) comprising adaptive antenna devices (AE) for receiving 
55 received signals (es) via time-variant channels, in which method 

• the received signals (es) from at least two adaptive antenna devices (AE) are weighted with antenna coeffi- 
cients (a) and are combined to form received data (z), 
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• error magnitudes (e) are determined in a device for determining error magnitudes (e) from the comparison of 
received data (z) and reference data (y) modeling these received data (z), 

• channel coefficients (RW) being provided for modeling the reference data (y), 

• these error magnitudes (e) are supplied to a channel estimator (ACE), 

• the channel estimator (ACE) jointly determines channel coefficients (RW) and antenna coefficients (a) mini- 
mizing the error magnitudes (e), and 

• the channel coefficients (RW) are provided for a detection of the received data (z) to be evaluated." (The 
detector is a Viterbi estimator whose short-time decisions are used outside the training sequence of a GSM 
cellular system data block). 

[0006] The RADEMACHER's LMS parameter optimization algorithm has recourse to three separate gradients in 
vectorial form and the optimization is carried out with respect to each gradient of the objective function, one independ- 
ently from the others, by using the same error (probably this originates the term "joint"); consequently inter-relationships 
between spatial and temporal terms are neglected and the opportunity to exploit these inter-relationships to improve 
the equalization definitively lost. 

[0007] A best profitable algorithm for joint estimating space-time radiomobile channel is disclosed in the International 
patent application, publication number WO 99/65160, inventor Umberto SPAGNOLINI et al, in the name of the same 
Applicant Siemens. The invention concerns an equalization method for the cancellation of isofrequential interferers 
coming from adjacent cells. The described algorithm exploits cross-correlation among space and time parameters for 
building the matrices used to perform the minimum error in the optimum criterion. 

[0008] The two mentioned inventions have been conceived for the channel estimate of a second generation cellular 
system and seem not immediately transferable to a third generation one to which the invention in subject is prevalently 
but not exclusively directed. Although the invention subject of the present application make use of joint space-time 
channel estimate, the solved technical problem that will be cleared later on is quite different from the mere joint space- 
time estimation of the two cited Application. 

[0009] Nevertheless from the above prior art the advantageous opportunity emerges to consider both spatial and 
temporal dimensions of the radio channel made possible by the exploitation of directive arrays from the BTSs. 
[0010] A modem design of a BTS receiver must unavoidable have recourse to sophisticated filtering and detecting 
algorithms, parallely to the prototyping phase of the receiver even a test of the joint operation of those algorithms 
becomes essential. This is generally performed by a computer simulation phase which needs a realistic modeling of 
the radiomobile channel. The modeling should take into account for Rayleigh flat and selective fading, shadow fading, 
path loss on Line Of Sight (LOS), Doppler effect, etc. Traditional channel models are described in the Recommendation 
GSM 5.05 issued from European Telecommunication Standard Institute (ETSI). These models are nevertheless only 
temporal and represent the channel by means of a certain number of time discretized pulse responses. A most recent 
model which best supports beamforming algorithms are referred to the COST 259 Directional Channel Model. This is 
a stochastic channel model for Space-Time (S-T) systems that includes both azimuthal and temporal dispersions. 
[0011] With reference to the Directional Channel Model depicted in fig.1 the channel response depends on N inde- 
pendent local scatters (clusters), each described by multiphats with associated amplitude, delay and azimuth. Ail the 
parameters in the model of fig.1 are independent random variables; the probability density functions are assigned 
according to the propagation environment. In addition to the line-of-sight cluster there are some other clusters, their 
numbers are Poisson distributed. 

[0012] The following four COST-259 radio environments have been proposed: Generalized Typical Urban (GTU), 
Generalized Bad Urban (GBU), Generalized Rural Area (GRA), and Generalized Hill Terrain (GHT). GTU has high 
time dispersion and low angular dispersion. GBU has high both the dispersions. GRA has low both the dispersions. 
GHT has low the time dispersion and high the angular dispersion. 

[0013] Considering an array of antennae having a steering vector a (9) (being 6 the incidence angle between a planar 
wave and the direction of the ULA), under the hypothesis of only one mobile station transmitting the channel pulse 
response estimable from a directive array is: 



where: c indicates the actual cluster; A/ C the number of clusters; N r the number of rays which the cluster reflects towards 
the BTS; 6 C , is a relevant direction and x C2 is the daiay of ray r C2 , and g (r) is the transmission pulse filtered by the 
reception filter. 

[0014] Since now the main problems concerning radiomobile channel estimation and uplink signal detection have 
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been duly considered as far as concerns the traditional second generation mobile systems, in case, equipped with 
directive arrays. That constitutes background notions for the invention that will be disclosed which is mainly directed 
(but not exclusively) to a third generation mobile system (UMTS). The usefulness of the discussed arguments depends 
on the fact that the new UMTS evolves from the preceding knowledge. However a simple translation is not enough 

5 because of the greater complexity of the UMTS in comparison with GSM due to the different multiplexing approach, 
complexity that is dramatically increased when an attempt to use directive arrays is considered. In the new UMTS 
ambit, differently from the GSM, more users are simultaneously active on the same carrier, and the BTS receiver shall 
provide for all of them by simultaneously receiving from all the antennas, so it will be soon realized what complexity 
comes out . In this case adequate innovations have to be introduced for the purpose of managing the complexity in a 

10 proper way suitable for optimizing the performances of the receiver while maintaining acceptable costs. 

[0015] Radiomobile systems based on Code Division Multiple Access (CDMA) are subjected to an International 
standard charged to a 3GPP study group. Presently there are several outstanding system proposals which differ each 
other mainly for the band width of the modulated carrier, the type of adopted full-duplexing either FDD (Frequency 
Division Duplexing) or TDD (Time Division Duplexing), by the use or less of an additional Time Division Multiple Access 

15 (TDM A), and finally for the synchronization of the BTSs. CDMA technique is implemented convolving each oversampled 
bit of the transmission symbols by a pseudo-noise spread sequence, temporarily assigned by the network both to a 
so-called Resource Unit (i.e. the Mobile Station) and the BTS for communicating full-duplex. The spread sequences 

being each other orthogonal, that means ideally can be seen as noise at the output of a matched filter to a specific \ 
sequence. The original band of the symbol sequence is widened consequently and the carried information scattered : 

20 on the whole transmission band. Midambles and data scrambling codes are used for building the transmission burst 

and both of them have one to one relation with the cell; they are semipermanent data for allowing the discrimination \ 
of the Resource Units (RU) of the serving cell from the interferents of the neighbour cells. As far as concerns the 
midamble, the KCDMA RUs must use K different midambles in order to be recognizable to the BTS, supposing the 
burst preliminarily synchronized. The Kmidambles are obtainable through Kphase shift of a unique periodic sequence 

25 having period P = KW chips, being Wthe maximum time dispersion of the channel. Obviously the greater is number 

Kthe longer the midamble will be. : 

OUTLINED TECHNICAL PROBLEM 

30 [001 6] The aim is that to outline the technical problem arising from the high parametrization of the channel estimation 
charged to a third generation BTS employing a directive array of M antennas to Space-Time (S-T) exploring the various 
channels. 

[0017] With reference to the Figures 2.3 and 4 the equivalent low pass discrete-time mathematical model of the 
uplink hybrid TDMA-CDMA mobile radio system is considered. In fig.2 an uplink transmission burst is depicted in which 

35 a Training sequence in midamble position is visible between two Data blocks 1 and 2, the last being followed by a l - 
Guard Period GP to avoid interference with the successive burst. The discrete-time model is obtained by sampling at 
the chip rate MT C the signals received by the antennas after the chip matched filter: Within the same cell, in the same 
frequency band and in the same time slot, K users are assumed to be simultaneously active. Each mobile station (MS) 
is equipped with a single transmitter antenna, at the base station (BTS) an antenna array of M elements is employed. 

^0 The users transmit in bursts each consisting of the two data blocks and the user specific training sequence (midamble) j 
for channel estimation. Each data block contains N s QPSK symbols of duration T s = QT C spread by a user-specific j 
signature c k of O chips for the fc-th user. The midamble consists of N m - N- W-\ chips, Wis the channel length. For 
example, in a first case referred to the TDD-UTRA specification N = K x W= 456 chips (K- 8 and W= 57); in a second 
case referred to a Chinese CWTS standard A/= 128 chips (K- 8 and W- 16). At the BTS the received signals are 

45 separated into two subsets of samples, the first one depending on the midamble training sequences and the second 
one depending on the two data blocks of symbols. The conventional approach for the estimation of the S-T channels 
from the samples of midamble is based on the Least Square (LS) estimate of the coefficients of the FIR filter that 
describes the channel model. Then S-T multi-user detection is carried out on the two blocks of signals received before 
and after the midamble in order to estimate the 2A/ s Kdata symbols in the burst. The detection needs the preceding 

50 knowledge of the channel estimates and the signature sequences for k = 1 K. It is known art that the midamble 

tails needs to be cancelled before data detection. ] 
[0018] The multi-user detection for space-time channels is known, for example, from the article of X. Wang, H. V. 
Poor, "Space-time multi-user detection in multiple CDMA channels", published in the review IEEE Trans. Signal 
Processing, Vol. 47, pp. 2356-2374, September 1999. This type of detection is useful to contrast with the so-called 

55 Multiple Access Interference (MAI). MAI interference depends on the fact that the K signatures are not perfectly or- 
thogonal at the receiver inputs regardless of the spreading codes used from the transmitter. The cause of MAI being 
both the dispersion of the delays on the MS-BTS channel (selective Rayleigh fading) and the imperfect chip synchro- 
nization. Preliminary to the multi-user detection a multi-user estimation shall be also performed; in this case channel 
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is estimated simultaneously for all the users so that ail the user are considered sources of signal for each other. 
[0019] Interceil interference further challenges the detection, specially when traffic increases and the cluster size 
reduction appears to be profitable for promoting a more intense reuse of the midambles and the scrambling codes. 
Obviously the multi-user detection is only charged to the BTS receiver because, differently from a mobile station re- 
s ceiver, it simultaneously manages signals from all the active users in the served cell, or in the sector in case of tri- 
settorial cells. 

[0020] Fig.3 preliminarily shows the various vectors and matrix of the system used for estimating the temporal channel 
M™) of the K active users receiving data from a single antenna m. 

[0021] Fig.4 extends the preceding system to the M antennas for estimating the complete S-T channel H = [H, 

10 H^) of the same K users. The meaning of the various elements involved, so as the dimensions of the vectors and 
matrices are clearly indicated in the Figures. 

[0022] With reference to the Figures 3 and 4 the propagation channel for the fc-th user is described by the S-T matrix 
that is composed of the M temporal FIR filters {h kn j : 

H, = [h M> ...,h kM ] r . (2) 

Each FIR h^ m represents the channel impulse response (of lenght W) for the link between the fr-th user and the m-th 
antenna. Here we assume M < W. The estimation of the S-T channel matrices {H k } is based on the transmission 
20 of the Ktraining sequences {X^} , where X k is the W x A/Toeplitz matrix that denofes the convolution with the /c-th 
midamble. Let Y =[y(1 ),..., y(A/)] b*e\he M x N matrix containing N samples of the signals received by the M antennas, 

H = [Hj hy and X = [X^ X rr are multi-user matrices suitably defined to include, respectively, the K channel 

matrices and the K training sequence matrices. The received signal can be written as 

25 

K 

y = y, H * x * - N = HX " N ' (3) 

30 where N = [n(1),...,n(A/)] models both the additive ambient noise and the inter-cell interference, it is assumed to be 
temporally uncorrelated and spatially correlated with space covariance matrix R n . 

The conventional approach for the estimation of the matrices {H^} is the direct LS estimation of the K x M x W 
coefficients of the FIR model (2). The residuals from the channel estimation are then used to estimate R„: 

35 a A A 

H=[H1 H K ] = YX H (XX") , (4) 

R„ = YY H - YX H (XX H )* 1 XY H . (5) 

40 

where denotes the Hermitian transposition. 

[0023] This approach does not use any a priori information about the S-T features of the propagation but it simply 
estimates brute force all the M FIR filters. This involves the following drawbacks: 

45 1 . a noisy estimate of the S-T channel in presence of the usual training sequences having necessarily limited length; 

2. a great complexity of the multi-user data detection because of the recurs to a full-dimensional space-time channel 
correlation matrix for all the K users. 

[0024] The second drawback is easily understandable, while as far as concerns the first one a relevant hint is done 
50 in the article of Monica NICOLI and Umberto SPAGNOLINI (also designated inventors in the present application), titled: 
"MULTIUSER SPACE-TIME CHANNEL ESTIMATION FOR CDMA UNDER REDUCED-RANK CONSTRAINT), pub- 
lished on IEEE 2000, page 147 - 151. The article has essentially theoretical nature and refers to the opportunity to 
constrain each matrix HW describing the KS-T channels to have Reduced-Rank (RR), in order to describe parsimo- 
niously (i.e. with minimum distortion) the S-T channel H = [H,,...^^. The rank is said to reflect the minimum number 
55 of orthogonal space and time channel that can be used under the parsimony constraint. The rank is further said to 
depend on the number of wavefronts that sampled by the array can be considered as independent. Therefore in a 
multipath propagation the rank-order depends on the angle of arrival and delay spreads compared to the resolution of 
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the array and of the signature waveforms (or their bandwidth). In a RR approach the MAI interference Is implicitly taken 
into account by estimating the Kchannejs (HW) k jointly. Further in the article (expression 9 pag. 147) the opportunity 
of spatially pr£-whitening the full-rank channel estimate has been highlighted in order to obtain the re- 

duced-rank {H W estimate. The performance of parsimonious channel estimation have been obtained and 

5 graphically presefffed*5y simulating the uplink of TDD-UTRA proposal for IMT-2000. 

[0025] The invention which is the subject of the present application fulfills a further step in the direction open by the 
article, namely that of really obtaining a space-time multi-user receiver for uplink, able to adaptively modulate the grade 
of parsimony in dependence on the variable channel condition. This is not so immediate from the preceding theory; 
besides the intrinsic complexity of the matter forces to solve further remarkable theoretical-practical problems faced 

10 to the aim of conjugate complexity with the true opportunities offered by the actual technology, taking the cost into 
account. All that will better appear in the course of the description. 

[0026] The following arguments further enlarge the parsimony criterion announced in the cited article. The accuracy 
to which a channel estimate can be performed using a known training sequence inversely depends (approximately) 
on the ratio between the number L unknown o1 parameters to be estimated and the length L trBining o\ the training sequence, 

is being said parameters, for example, taps of the temporal and spatial filters. Recurring to training sequences very length 
it should be possible estimate the L unknown parameters near exactly, because the variance of the estimation decreases. 
Choosing to estimate a lower number of parameters (L' unkn own < L unknowri ^ e estimate should be run into an error, 
also termed polarization, due to the under-parametrization of the channel. Nevertheless in this event the estimation of 
the remaining parameters should be better, because the ratio L' un knowr/ draining iS lower than L U nknowr/ ^training- ^ e * u " 

20 rank estimation of the channel matrix H requires K x M x W parameters, while a Reduced Rank estimation requires 
a lower number depending on the skill in the exploitation of the a-priori information about the space-time propagation 
features along the radio path. Therefore in presence of Ltmimng finite, a trade-off between under-parametrization and 
accuracy of the estimate shall be pursued. A consequent parsimony principle can be issued to state that "it doesn't 
need estimate more parameters than the strict necessary". 

25 [0027] Fig. 5 graphically shows the parsimony principle relatively to a hypothetical RF scenario. The abscissa repro- 
duces a scale proportional to the parameters to be estimate and the ordinate the MSE error of the expression ||H - A H|P 
/ ||H|| 2 (being A H the estimated channel) which represents the accuracy of the estimate. In the hypothetical RF scenario 
on which the Fig.5 is based, the minimum MSE is obtained in correspondence of a parsimonious estimation using 3/8 
of the maximum parameters. The MSE on the upper curve is the sum of two terms: the distortion error, that is large 

30 for low rank and depends on the characteristics of the propagation channel (delay/angle spreads, number of clusters), 
and the noise that increases with rank and depends on the interference level, other than the increasing of L unknowr / 
^training- 

[0028] The application of the parsimony principle in the design of a multi-user receiver for uplink TD-CDMA involves 
some practical implications that could at first-sight discourage the design of the receiver; in fact, the parsimonious 

35 estimation of the channel seems not so parsimonious like the word indicates, that because the reduction of the number 
of spatial and temporal filters needs in any case the preliminary full-rank estimation of the multi-user channel matrix 
and the additional handling for the rank reduction. Nonetheless in spite of a slight increment in the computational 
complexity of the channel estimation, it has been considered above that the rank-reduction allows a considerable 
performance as far as concerns the MSE error of the estimation. The complexity of the multi-user detection is instead 

to always reduced in the reduced-rank detection, because the number of involved spatial and temporal filters is lower. 
[0029] Turning the attention to the two firstly cited inventions concerning GSM, it can be argued they have not a so 
intense parametrization to force solutions based on parsimonious estimation, although parsimony should be even 
exploited in a GSM-BTS. As far as concerns the article of the M. NICOLI and U. SPAGNOLINI, some relevant sugges- 
tions are instead taken, but relevant new features have been added to solve a technical problem apparently not yet 
well outlined in the new UMTS context. The reduced-rank estimation method disclosed in the article appears "static" 
in the choice of the optimal rank and could be further improved to allow optimal estimation. That because further 
investigations (that will be confirmed discussing the Figures 10 and 11 about the variation of the MSE error (Fig.5) for 
the various ranks in function of SNR (Signal to Noise Ratio) or SIR (Signal to Interference Ratio), have proved that: 
for low SNR (or SIR) low rank is the preferred solution as it has the least number of unknowns to be estimated. For 

50 large SNR the distortion (polarization of the estimate) becomes remarkable and more ranks are needed. The tracking 
of the best structure of the S-T channel matrix seems then a difficult task to be performed by the receiver, mainly 
because said structure shall cope with the multipath structure of the radio propagation depicted in the expression (1). 

PURPOSES OF THE INVENTION 

55 

[0030] The main purpose of the present invention is that to indicate a multi-user receiver, preferably for uplink CDMA, 
performing joint space-time multi-channel estimation and joint multi-user data detection, able to adaptively estimate 
the optimum trade-off between distortion and noise variance as dependent on the interference level and the propagation 
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environment. That also addresses the purpose of the invention to indicate a modular S-T receiver whose structure is 
upgraded or downgraded consequently for estimating the actual leading space-time characteristics of the channel. 

SUMMARY AND ADVANTAGES OF THE INVENTION 

[0031 ] To achieve said purposes the subject of the present invention is a Base Station receiver, as disclosed in claim 1 . 
[0032] From the above claim it can be appreciated that the receiver is usable in all the cases where the S-T channel 
can be preliminarily estimated by any method (e.g. LS channel estimate) and a re-estimation procedure can be applied 
to reduce the number of unknowns to only those that are effectively necessaries when compared with channel and 
interference conditions. The advantages of the invention are evident in TD/CDMA systems for multi-user approach, 
where the complexity is high, but it is applicable also in the GSM system for single user approach. 
[0033] In comparison with the conventional approach where the channel model is assumed to be full-dimensional 
(Full-Rank), here the complexity of the model is reduced and the number of parameters to be re-estimated for the 
single user becomes r^M - W), usually « MW. The approximation with realistic propagation channel model COST- 
259 for the TDD-UTRA standard holds as far as r k ranges from 1 + 4 and the array is composed of M = 8 antennas. 
[0034] It should be noticed that at least one space and time filter are used and thus r k > 1 , this implies that the 
estimation structure is modular in the sense that the number of cascaded space and time filters can be increased by 
starting from one space-time filter up to M space-time filters. The actual number r k \s adaptively estimated by optimizing 
the trade-off between distortion and noise variance as dependent on the interference level and the propagation envi- 
ronment. This reduction of the dimensions of the problem results not only in better estimation performance when eval- 
uated in term of mean square error (MSE) but also in the reduction of the complexity of the space-time detector. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0035] Further objects and advantages of the present invention will be made clear by the following detailed description 
of an embodiment thereof and the annexed drawings given for purely non-limiting explanatory purposes and wherein: 

• Fig.1 shows the channel model COST 259. 

• Fig.2 shows a typical transmission burst for TD-CDMA 3GPP radiomobile systems. 

• Fig.3 shows the basic vectors used for estimating multiple temporal channel with a single antenna. 

• Flg.4 shows the basic matrices used for estimating a spatial-temporal channel with M antennas. 

• Fig.5 shows a graphical representation of a principle which inspires the channel estimation of the invention. 

• Fig.6 to 9 show the architecture of the Base Station receiver of the present invention; 

• Figures 10 and 11 show families of curves reproducing comparative simulation results of the channel estimator 
of fig.6 

• Figures 12 and 13 show families of curves reproducing comparative simulation results of the multi-dimensional 
BASE STATION RECEIVER of fig.6. 

• Fig. 14 shows families of curves reproducing comparative simulation results of a mono-dimensional BASE STA- 
TION RECEIVER of fig.6. 

DETAILED DESCRIPTION OF AN EMBODIMENT OF THE INVENTION 

[0036] Figures 1 to 5 have already been discussed. Before the description of the preferred embodiment the following 
theoretical arguments are proposed for justifying a successive so-called pre-whitening of the estimated matrix (4), in 
order to make the K users each other spatially uncorrelated. For this aim it is needed to evaluate the covariance matrix 
of the channel estimated as: 



R„ h = El5h5h H ] = (XX H )- 1 <8>R n (6) 

A A 

where Sh =vec(8H) = vec(H - H) denotes the stack of the columns of the estimation error matrix H - H into a vector and 
<8> denotes the Kronecker product. Some simplification are in order only for the reduced complexity estimation receiver 
that will be short-term disclosed: 

1 . The covariance matrix of the interference can be approximated from the estimate (5): 
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(7) 



2. The inverse of the training sequences for different users is computed by assuming that the training sequence 
are approximately orthogonal: 



10 



XX 



H _ 
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x x xf x*x? X A -X£ . 



x,xf 0 
0 x 2 xf 
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0 



Xk-X/J 



(8) 



and therefore, for the temporal pre-whitening, it follows: 
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(9) 



3. The covariance matrix can be approximated as 



25 
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0 (X 2 X.f )" ! « Rn 



0 
0 



(Xa-X^)- 1 ®^ _ 



(10) 



and therefore for the pre-whitening the estimates are decoupled thus yieldind 



El5h5h£] - (X k X?)" 1 ® R„ (11) 

[0037] Fig.6 shows the high level architecture of the Base Station Receiver of the present invention. With reference 
to the figure the receiver is constituted by the following blocks numbered from 1 to 12: ULA, Front-end, Burst-Mem, 
40 Training-store, Channel model, Channel estimator, R„ model, Data Receiver, which in its turn includes: Decorrelator, 
SWF & S-T Matched filter, and spreading code. The names of the blocks indicate their meaning. M-antenna Uniform 
Linear Array 1 is connected to the Front-end 2 which provides to: 



• correctly receive the RF signals transmitted from all the K users, 

• filtering and down-converting the RF signals, 

• Q-PSK demodulating, sampling and A/D converting the base band signal, 

• base band filtering with a digital filter matched to the chip obtaining the base band digital samples Y corresponding 
to the sum of the K transmitted bursts (see 3) that are stored in the block Burst-Mem 3 to be processed by the 
base band receiver of the invention. 

[0038] The receiver processes the composite signal Y in two alternate steps ideally set by a switch 4: a first one 
involves blocks Training store 5 and Channel estimator 6, and provides to estimate the channel parameters stored in 
block Channel model 7 and R n model 8, while the second one uses the preceding estimates for detecting the two semi- 
bursts of data parallely for the K users which keep the Base Station Receiver busy. Block Channel estimator 6 receives 
the samples Y corresponding to the transmitted training sequences and the sample X correspond ing^to the stored 
training sequences and estimates both the channel matrix H k and the spatial noise covariance matrix H n , which rep- 
resents the content of block Channel model and R n model respectively. The two estimated matrices are sent to the 
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block Data Receiver 9 referenced in Fig.8. 

[0039] Fig.7 shows the detail of block Channel estimator 6 of the preceding figure. With reference to the Fig.7 the 
outlined blocks are listed from 13 to 22, they are the following: LS-ST channel estimator, Spatial covariance matrix 
estimator, Temporal whitening, Spatial whitening, REDCOMP, selector, VAR-DIM, RECOLOURING, R n re-estimator, 
and Memory. Some blocks and signal are indicated with subscript k to be referred to the user k* h , in such a case it is 
obvious that K similar blocks and signals are provided even if not depicted for reason of brevity. 
[0040] Block LS-ST channel estimator 13 is reached from the signals X and Y (fig.6) and exploits them to firstly 
perform a full-dimensional Least Square (LS) estimate H^of the S-T channel in a known way, for example by means 
qf Fast Fourier Transform (FFT) and its Inverse. Block Spatial covariance matrix estimator 1 4 performs the estimation 
R„ of a spatial covariance matrix R„ starting from the received signals X and Y (fig.6). The estimated channel is 
temporal whitened by the Temporal whitening filter 15 and then forwarded to the spatial whitening filter 16. For the aim 
of temporal whitening signal X k relative to the user feeds the input of Temporal whitening filter 15 by means of a S/ 
P converter (not shown in fig.7) which parallelizes upon Wways the single chip of the sequence relative to the k* h user. 
[0041] The channel model which comes out from the complex chain of algebraic operators is a spatial and temporal 
whitened channel model named H'^. This model is full-dimensional and needs to be complexity reduced for the reasons 
outlined in the introduction. Blocks REDCOMP 17 and VAR-DIM 18 together provide to this aim. Processor block 
REDCQMP is the most relevant of the two, it fulfills two tasks: a first one is that of decomposing the whitened channel 

model of the ffl 1 user into a set of filtering cascades A^B'^. each constituted of a beamformer A^ = [a'^ ,a'* (ri J 

followed by a temporal FIR B* k = [W^ b\ f J of dimensions Mx r k and Wx respectively. Each space and time 

fitter are constrained to be orthogonal (after pre-whitening) by the particular operation of the REDCOMP block which 
will be explained short time: 

i u* n = n f or m ± n 
I b£A.n = 0 

A, 

The criterion performed in the decomposition indirectly seeks for the minimum error between channel model and 
the new representation A H B H which is better idoneous for complexity reduction. Complying with the decomposition 
criterion allows the re-estimated parameters of the various filters be determined, but the number of the filtering cascades 
other than a first one is still pending. So the second task charged to REDCOMP processor 17 is to determine the 
effective number r k « M of the filtering tandem relative to the A" 1 user, in accordance with the above mentioned par- 
simony principle limiting at the only effectively needed the channel parameters to be re-estimated. The complicated 
operation of processor REDCOMP is supported by the memory block VAR-DIM 1 8 which stores the variable configu- 
ration of the reduced-complexity model during its elaboration before the exact number r k be determined. Inside the 
block VAR-DIM a number r^of filtering cascades are visible, in which: beamformers relevant to the spatial components 
a . are denoted by s h while temporal FIRs relevant to the temporal components b . are denoted by t t Only the first 
filtering branch is drawn in solid, because the relative presence is always assured, &e remaining branches are dotted 
because their exact number depend on the result of the second task charged to the processor REDCOMP 17, in its 
turn depending on the actual level of the inter-ceil interferences. For this aim block "selector" 19 selects the branch of 
block VAR-DIM 18 which has to be presently built up; after which a series of arrows superimposed to the various FIRs 
select, in turn, filter s,- and f/Of the selected branch whose components a and b . have to be updated and stored. 
Vectors that result from the equivalent branch filtering are summed up for'all the r£ branches as in the following: 

£>;X" = a^"-h;. (13) 

and the product (13) is both delivered to the processor REDCOMP 17 a^id to the block RECOLOURING 20. The first 
one performs the error computation during re-estimation of the channel while the second one performs spatial and 
temporal inverse whitening filtering to obtain the recotoured product A^B^, which is stored in the Memory 21 at the 
disposal of the Data Receiver of Fig.8. Lastly, the Channel estimator in Figure shall re-estimate the spatial noise cov- 
ariance matrix R m directly in the form fc^ suitable for spatial whitening, making use of the new reduced complexity 
channel parameters a^and b k h This task is performed from the block R n re-estimator 22 which receives the recoloured 
product A^B H and signals Y and X fc and returns . Also the re-estimated covariance ma\nxK^ H/2 is stored in memory 
21 for the successive use in Joint Detection. Block memory 21 shows, for simplicity, a dashed box representing room 
for the spreading sequence C* of fig.6 assigned to the K users. 

[0042] Now some mathematical hints are done to get an insight into the Channel estimator of Fig.7. The opportunity 
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A r k X k 

of approximating the Full Rajjk H k with r k spatial signatures {a k J and r k temporal filters {b k J , in such a way that 
redundancy of the Full Rank H k '\s avoided in the latter representation, is achieved by constrainingftrlat each beamformer 
and temporal FIR filter is orthogonal to the other r k - 1 filters. Orthogonality is constrained for pre-whitened estimates 
in order to take into accounts the spatial correlation of the interference and the non-ideality of the training sequences 
(i.e., XX H * I). In comparison with the conventional approach where the channel is assumed to be full-dimensional (r k 
= M, A k = \ M and B*= H k ) t here the complexity of the model is reduced and the number of parameters to be re-estimated 
becomes r^M + W), usually « MW. It should be noticed that the number of degrees of freedom remains the same as 
the constraints (1 2) corresponds to - 1 ) equations and the unit-norm constraint are r k equations thus leading to r k 
(M+ W) - r = Oc(M+ W- r k ) parameters to be re-estimated. 

[0043] Since the users can be considered as decoupled, for the /c-th user the channel estimate is whitened as 

H*=R-„ H/2 H,(X^) H/2 (14) 

this is approximated by the low-dimensional channel A B H , where A = [di k V ...,a krk ] = R A^ are r k orthonormal 
whitened beamformers and B^ = [b^ ^ , .-.b^ ] = (X k X^v /2 B^ are r k orthogonal whitened FIR filters: 

R*«A fc B'? (15) 

The ^-dimensional spatio-temporal components (Reduced Rank) is obtained by minimizing the error of the approxi- 
mation: 

{A; ; B'J = arg min \\R k - A k B' k H \\ F . (16) 

A:.,B' b 

in the Frobenious norm. The optimization is known to be given starting from defining the following matrixes: 

A A, A 

Z k = H k H, H , (17a) 

for n = 1 , ... r k 

0 k =eigv rk {Z k } (18) 

A 

where eigv^ denotes the M x r k matrix of the r k leading eigenvectors of the matrix Z*. Whitened estimated A fc , B^ are 
given by: 

A\ = 0, (19) 
B'k = (X k X^) 1/2 H^R n 1/2 0 k (20) 



A definitive estimate of ST channel is given by: 



A*=CA* (21) 



BfcMX^tVX (22) 

With respect to the conventional estimation this approach provides an improvement in the estimation performance as 
well as a reduction of the complexity for the equalization processing, some relevant results are shown in Figure 1 0 to 
14. The space noise covariance matrix for the data estimation is re-estimated from the residuals of the channel esti- 
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mation N = Y - Z. K ^ A^B H p(^ as: it, = NN H IN . 

[0044] The number r k of space-time filters inside VAR-DIM block 18 is not known and need to be estimated by 
increasing approximations of the reduced-complexity re-estimate. It should be noticed that at least one space and time 
filter are used and thus r k > 1 , this implies that the estimation structure of block VAR-DIM (fig.7) is modular in the sense 

s that the number of space and time filters can be increased by starting from one space-time filter up to M space-time 
filters, according to the particular operation of the Reduction complexity processor REDCOMP 17. The optimum r k is 
the one that minimizes the mean square error (MSE) of the channel estimate \\H k - A A B H || F but this MMSE optimization 
cannot be carried out in practice. Here r k is obtained according to the information theory criterion based on the Kullback- 
Leibler information measure as proposed in M. Wax. I. Ziskind, "Detection of the number of coherent signals by the 

10 MDL principle", IEEE Trans. Acoustic Speech and Signal Processing, Vol. 37, pp. 1190-1196, August 1989, which is 
incorporated by reference, for the estimation of the number of independent sources impinging on a uniform linear array. 
The information theory criterion is based on the Minimum Description Length principle (MQ|_) , applied in the present 
case for the description of the temporal covariance matrix of the full-dimensional channel H relative to the k? h user. 
The MDL algorithm for the selection of the optimum r k inside the modular structure of block VArt-DIM 1 8 is implemented 

is by REDCOMP processor 1 7 through the following steps: 



1 . define a vector L(M x 1 ) of null element; 

2. compute the matrix T = H^H H and its eigenvalue (/= 1, M)\ 

3. execute for p = 0 : M - 1 
20 4. 



25 



L(p-1; = -W(M—p)log 



' (n^)*'" 



-p(2M-p-l)\o*(W) (23) 



5. minimum of vector L -> (pos. vat) (position in the vector, minimum value) 
30 6. if (pos- 1) = 0 then r fc =1 else r k - pos -1 

7. r k is the output of the MDL algorithm (repeated for each user) 

It should be noticed that r k > 1 (at least one ST element is adopted by the receiver) and the final value of r k 
corresponds to the minimum MSE on the upper curve of fig.5. 

35 [0045] At this point of the description the reduced-dimensional (Reduced-Rank) channel model A^B^ recoloured, 
and the re-estimated space noise covariance matrixa^ w/2 are stored in memory 21 , so that data detection (estimation) 
performed by block Data Receiver (fig.6) starts by turning the switch 4 (fig.6) in the relative direction. 
[0046] With reference to the fig.8 it can be appreciated that the Data Receiver 9 (fig.6) includes the cascade of the 
following three blocks: SWF 11 a, as many parallel S-T matched filters 11 b as the K users, the DECORRELATOR 10 

*o in the manner of a multi-user detector. The memory block 21 (fig.7) is dotted in fig.8 with enabled room (in solid) for 
the signatures C^for all the K users. In this processing step of the memorized burst the transmission of one single data 
block (before or after the midamble) is considered. The processing step avails of the preceding estimated elements: 
R^ w/2 and re-estimated channel A^B*. An additional block 23 named RT-S-WT is cascaded to the only A k beamformers 
stored in memory 21 for performing a spatial whitening of A k by using the re-estimated covariance matrix ft,. More 

*5 precisely this whitening is obtained through the following espression: 



A* = R^ /2 A fc (24) 

50 

Equation (24) contains all the at.i,.... a*.,. ! beamformers. 

[0047] Following is described the discrete-time model for the detection of a data transmitted by the K users. Let d k 
(/) denote the f h symbol transmitted by the fc^ user, d(/) = [d,(/> — c^/)] 7 is the vector containing the I th symbol for all 
the K users, d = [d(1) r - d(A/ s ) 7 ] r is the overall data vector of length N S K. For user k the convolution of the spreading 
55 sequence c k with the € th temporal channel b k € is denoted by the composite signature s k ( = C^b^ c of length Q + W - 
1 , with € = 1 , .... r. For simplicity of notation all the users are assumed to have the same number of S-T dimensions r k 
- r, V/c. The multi-user Data receiver estimates the user data d from the whole received signals y(n) = [y t (/) -y^/)] 7 , 
with n = 1 N S Q - W- 1 , using the beamspreading filters^ and the signatures s^. 



11 



EP 1 255 387 A1 



[0048] From the discrete-time model of above, the tf h S-T matched filter block 11b can be qualified as including r 
beamformers 



10 



20 



25 



for € = 1 ,...,/* (whitening space matche^ filters). Each beamformer is followed by a corresponding temporal filter s ( e . 
The outputs of the rtemporal filters s^ f are sub-sampled by factor O before reaching the r inputs of a summation iblock 
which provides the filtered sequence WO 1 matched to the estimated channel. 

[0049] The complete output y(0 - fo(0 • • faWF -of the K space-time matched filters 1 1 b is obtained by whitening the 
received signal y(n) with the spatial filter SWF 11a which executes the productft^Y, then the whitened signal is 

passed through the Kr beamformers «£, = {f£ m *k*) H ; for k- 1 and t = 1 r(whitening space matched filters). 

Each beamformer a£< is followed by a corresponding FIR temporal filter^. The filtered sequence j/*(0 is derived by 
summing up the outputs of the s H filters for i = 1,.., r and the filtered outputs for all the K user- 
sy(i; = ■ • ■ {//<f0l T are fec^into the DECOR RELATOR block 10 which acts like a linear S-T multi-user de- 

tector that estimates the user data d(/) = [</,(/) — c/^/)] 7 . 

[0050] It can been observed that this receiver has a simpler structure with respect to the full dimensional space-time 
receiver, as it does not exploit all the spatio-temporal dimensions in the matched filtering. Indeed for the conventional 
receiver the structure is the same of Fig.8 but with M space-time matched filters instead of r, a ke is a M-dimensional 
vector of zeros with a one entry in the €-th position (A k = \ M ) and the temporal matched filter has b^ -h ke (B k =H^. 
[0051] Fig.9 shows the linear DECORRELATOR block 10 based on an advantageous embodiment founded on 
Cholesky factorization method. With reference to the Figure the proposed general architecture of the DECORRELATOR 
10 includes the following blocks: R s computation 24, R, computation 25, Hadamard multiplication 26, Cholesky factor- 
ization 27, equation solver A 28, and equation solver B 29. For clarity, the same memory block 21 of Figures 7 and 8, 
and block RT-S-WT 23 of fig.8 are reproduced. Mathematical arguments supporting the diagrammatical representation 
of fig.9 begin from arranging the matched filter outputs of fig.8 for all symbols into the following vec- 
tor^ - [y T U; •y r (- v J )! r , then the system model for the overall data block reduces to: 



30 



y == Rd -r n 



(25) 



35 



45 



where the noise vector is n -vVi(0, R). R is the overall S-T correlation matrix with dimensions (N S K x N $ K) and has 
the block Toeplitz structure denoted by 



R-Tv,(fR (0) R<*>]), 

' R(0) ' R (l) 



R = 



... R< A > 

Rl-U R(0) RlD ... R(A) 

* - \ 

R i-4) Rl -i) R ,o) R m 

R (-A) R (-0 R w 



(26) 



where A^ 2 ^,' is the temporal spread in symbol intervals. Each block RM = R<-'> H has dimensions Kx K(i= 0 A) 

and can be written as the Hadamard (element-wise) product between the spatial and the temporal correlation matrices: 



50 



oR!*G 



(27) 



55 



G = Ik®1 p 1 f is a r-dimensional vector of ones and o denotes the Hadamard product. 

The Kr x Kr matrix R s contains the correlations between the Kr spatial beamformers^,', whereas the Kr x Krmatrix 
R <y * depends on the correlations between the temporal filters s k( (delayed of /A). 

As* already said the outputs y of the S-T matched filter 11b (25) are fed into the DECORRELATOR block 25 that 
estimates the user data d(/) = [<*,(/)•■ •d fc (/)] 7 '. The ideal detector removes both Inter-Symbol Interference (ISI) and Multiple 
Access Interference (MAI) by processing the whole received data block a - M y. In particular, the decorrelating de- 
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tector is M = R- 1 , the MMSE detector is M = (R+a 2 l WK )* 1 if the data symbols are independent and uniformly distributed. 
[0052] The space-time structure of channel estimate allows to easily compute the ST correlation matrix R: the element 
by element multiplication that results in each block evaluation has low cost compared to conventional receiver 
(especially with r « M). By Cholesky decomposition of matrix R it is possible to easily resolve linear system that supply 
5 data estimation. The linear equation solver get out estimated data 3 that are input for successive channel decoding 
(not shown in the fig.9). 

[0053] Most of the computational complexity of the linear system resolution comes from the Cholesky factorization 
of the sparse and large correlation matrix. As R has a block Toeplitz structure, this factorization requires in the order 
of Ot/^NgA] operations as shown by J. Rissanen, "Algorithms for triangular decomposition of block Hankel and Toeplitz 

10 matrices with application to factoring positivew matrix polynomials", Math. Comp., Vol. 27, pp. 1 47-1 54, January 1 973. 
If the size of the data block N s is very large ideal detection may not be feasible. A more practical solution can be derived 
by using a sliding observation window that spans N a « N $ symbol lengths, as it results from M. J. Juntti, B. Aazhang, 
"Finite memory-length linear multi-user detection for asynchronous CDMA communications'*, IEEE Trans, on Comm., 
Vol.45, pp. 611-622, May 1997. An alternative solution is an approximation of the Cholesky factor that exploits the 

f 5 block Toeplitz structure of R and the property that RM = 0 for |/| > A, with A « Ng. It turns out that the Cholesky factor 
R 1/2 is nearly block Toeplitz and can be calculated from the factorization of a N a K x N a K matrix of the cited Rissanen. 
For both methods the computational complexity is reduced to ty^N^], with a negligible degradation in performance. 

SIMPLIFIED EMBODIMENT 

20 

[0054] Single space/time receiver is based on the choice r k = 1 (i.e., no MDL is performed). Below is the table (values 
are x 10 6 ) with computational complexity of standard receiver based on Full-Rank (FR) channel estimate (4) and the 
reduced complexity estimate with r k = 1 . The values are the number of real-multiplications (MUL) for each of the steps 
of the receiver: channel estimation, S-T matched filter, computation of the correlation matrix (26), Cholesky factorization 
25 of correlation matrix (26), forward/backward substitutions. The comparison is based on the following values: M= 8, K 
= 8, W= 57, Q = 16, N s = 61 (for each of the semi-burst). The Cholesky factorization is based on the choice N a = 7 
(the factorization for N a = N s = 61 would have implied 1.3 • 10 6 MUL). 





H est 


R„ est ] ST MF 


Corr. matrix Oecorrelation 


Total 


FR 


0.27 


■5.31 2.2 


0.4 | 0.62 


3.3 


r* = l 


0.37 


6.32 ! 0.81 


0.051 | 0.62 


3.2 



35 Since the eigenvector for the largest eigenvalue has to be computed according to the equations (1 7a-18), the eigen- 
vector computation is based on the power method (here 3 iterations are considered as enough for the algorithm to 
converge in many cases). Here the computation of the complexity is based on the approximation I which holds 

true for TDD-UTRA standard. 

[0055] As a further comment the conventional and the reduced complexity method have comparable computational 
40 complexity as the channel of conventional method is based on M FIR filters of W samples each is reduced to one FIR 
filter with W samples and one beamformer with M samples. This simplify the correlation matrix (26) in the multi-user 
detection. 

PERFORMANCES 

45 

[0056] In the remaining Figures the performance of low-dimensional channel estimation and S-T Multi-User Detection 
(MUD) are evaluated by simulating the uplink of the UTRA-TDD described in standard ETSI, "Universal Mobile Tele- 
communications System (UMTS): Physical channels and mapping of transport channels onto physical channels (TDD)" 
ETSI, TS-125 221 V3.3.0 2000-06. Numerical results are for M= 8 omnidirectional antennas, K = 8 users are active 

50 in the cell, the channel-length is W- 57. The traffic burst 1 of TDD-UTRA is considered, spreading codes have O = 
1 6, raised cosine pulse with roll-off 0.22 and chip rate 3.84 Mchips/s. Data blocks of N s = 61 QPSK modulated symbols 
are considered, training sequences are chosen according to standard specifications with length N m = 512. The prop- 
agation channel is simulated as defined by the Generalized Typical Urban environment (GTU) of COST-259 Directional 
Channel Model. In the simulations both spatially uncorrected and correlated noise (generated by 3 intercell interferers) 

55 are considered. 

[0057] Figures 10 and 11 respectively for uncorrected and correlated space noise, both compare the performance 
of the low-dimensional channel estimate with the full dimensional estimate in terms of the normalized MSE |]AH|| 2 /||H|| 2 
vs. the signal to noise ratio SNR= E[(||h(*' m )|| 2 )]/a 2 . The considered case are: Reduced-Dimensional channel (RD) with 
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fixed r, RD channel with estimated r, and Full Rank (FR) channel. As expected for low SNR the 1-D approximation is 
the preferred solution as it has the least number of unknowns to be estimated. For large SNR the distortion becomes 
remarkable and more space-time dimensions are needed. The channel estimate with MDL selection of the number of 
space-time filters (circle line) outperforms the fixed-dimensionality estimates (with r = 1,2,3,4, dashed lines) and the 
FR estimate (thick line) for all the SNR values. 

[0058] The receiver, complete with channel estimation and S-T MMSE multi-user detection, is considered in the 
remaining Figures. Fig.12 shows the performance for uncoded bits vs. SNR = E bi /N Qt for uncorrelated space noise. 
Fig. 13 shows the performance of the receiver in presence of 3 intercel! interferers evaluated in terms of BER vs. signal 
interference ratio (SIR). The legend is again the previous one. From the curves emerges that, if a fixed number of filters 
is adopted for the low-dimensional estimate (dashed lines) the distortion becomes very marked for large SIR/SNR, 
and the estimate performs worse than the FR (thick line); the adaptive selection of the dimensionality by MDL criterion 
(circle line) guarantees a gain with respect to FR of approx. 3dB for all SNR/SIR values. 

[0059] In Fig.1 4 the performances of a mono-dimensional receiver (r = 1 ) are evaluated for a propagation channel 
with lower spatio-temporal dispersion. The channel model is GTU, with number of cluster fixed and equal to 1, K ( = 3 
inter-cell interferers are simulated, the signal to noise ratio is SNR = 20dB. The 1 -D receiver provides a gain of approx. 
3dB even for large SIR. 

[0060] Simulations results presented in the Figures 1 0 to 1 4 have proven that the Base Station receiver of the present 
invention provides improvements in system performance (evaluated in terms of BER) by at least 3dB in signal to noise 
ratio (SNR). That because a particular channel estimation has been adopted by which the number of filters is selected 
adaptively according to the interference level and the characteristics of the propagation environment. 



Claims 

1. A Base Station receiver connected to a directional array of M antennas for receiving bursts of data from K mobile 
stations of a cellular telephony system, the receiver having: 

A 

• a full-dimensional space-time channel model (1 3) in form of memory matrix table (H) including a certain number 
of elements corresponding to relevant W samples of the dispersive length of the channel pulse response of 
all the K users, with k = 1.....K, for the link between the tf* mobile station transmitting antenna and the nt h 
receiving antenna of the Base Station, and for all the N samples of the tf* training sequence assigned to the 
tf" user, being the K training sequences each other approximately uncorrelated; 

• multi-channel parameters estimating means (6) for estimating the whole samples of said full-dimensional chan- 
nel model (1 3) for the K users simultaneously by correlating the base band digital samples of the signals 
received on the array with the corresponding K training sequences (X) of samples stored (5) in the receiver, 

• multi-channel noise and interference estimating means (1 4) receiving base band digital samples (Y) of the 
signals received on the array to which subtract the residual of the channel estimation for obtaining an estimate 
of the spatial noise for all the K users, being the temporal correlation of noise and interference considered as 
known; 

• means for convolving (1 1 ) base band digital samples of data transmitted from alt the K users with relevant 

estimated parameters of the channel model , obtaining K" sequences jf* *£) of the K users spatially 

and temporally matched to the channel; 

• multi-user detecting means (1 0) receiving the estimated parameters of the channel model, the K matched 
filtered sequences of the K users for jointly de-correlating the original data .... ) transmitted 
from all the K users, 

characterized in that said multi-channel parameters estimating means (6) further includes: 

• reduction-complexity means (1 7, 1 8, 1 9) for re-estimating said f ull-dimensiona space-time channel model ( 1 3) 
with a reduced set of re-estimated parameters (A^, B'") in order to obtain one or more cascades (s t , f 1t 

s r t r ), each cascade including a beamformer(s 1 s r ) and a temporal FIR filter (t, f p ), being said cascades 

decoupled each other and their numbers depending on the specific user; 

• reduction-complexity control means (1 7) arranged for limiting the re-estimation to a number r k ot said cascades 

(s v f t s r f r ) for each of the K users, r k being chosen dynamically as the number of independent cascades 

that minimizes the redundancy of the full-dimensional space-time channel model (13) according to the infor- 
mation theory criterium. 

2. A Base Station receiver according to claim 1 , characterized in that each of said cascade (s v ... s n t r ) is equiv- 
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alent to a user branch matrix, and for each user equal elements in the associated r k branch matrices are summed 
up to obtain complete cascade matrix; 

3. A Base Station receiver according to claim 1 or 2, characterized in that it further includes means for constraining 
5 temporal uncorrelation £1 5), ajso teamed temporal whitening filtering means, acting upon the estimated full-dimen- 

sional channel model H = [H„ ...,H k , ...,H K ] (13) before using it for the re-estimation, for obtaining a temporal 
whitened channel model (15) by uncorrelating said model (13) from the intrinsic temporal correlation residual from 
the non ideality of the training sequences (X k ). 

10 4. A Base Station receiver according to claim 3, characterized in that said means for cpnstraining temporal uncor- 
relation (15) being arranged for multiplicating said channel model estimated matrix by a first temporal factor 
(X A X H ) W2 where X k indicating the tf A training sequences (X^) of samples stored in the receiver (5) and said expo- 
nent vindicating transposed coniugated operator. 

15 5. A Base Station receiver according to one of the preceding claims, characterized in that said mufti-channel noise 
and interference estimating means (1 4) being arranged for estimating a space covariance matrix R n of the sample 
collection of additive ambient noise and inter-cell interference among the various antennas, and for storing it in a 
noise channel memory (21 ) of the receiver. 

20 6. A Base Station receiver according to claim 5, characterized in that it further includes means for constraining 
spatial uncorrelation (1 6), also termed spatial whitening filtering means, acting on said temporal whitened channel 
model (15) before using it for the re-estimation, for obtaining a temporal and spatial whitened channel model 
(1 6) for the k* h user by further uncorrelating said temporal whitened channel model (1 5) from the spatial correlation 
due to the inter-cell interferes. 

25 

7. A Base Station receiver according to claim 6, characterized in that said mear>s for constraining spatial uncorre- 
lation (16), being arranged for using the estimated space covariance matrix R n by multiplicating said temporal 
whitened channel model (1 5) for the tf A user by a first spatial factor R"*' 2 , where exponent H indicating transposed 
coniugated operator, and obtaining said temporal and spatial whitened channel model H (16). 

30 

8. A Base Station receiver according to one of the preceding claims, characterized in that said reduction-complexity 
control means (17, 18, 1 9) being arranged for re-estimating said temporal and spatial whitened channel model H'^ 
(16) for the k* h user for obtaining said reduced set of re-estimated parameters (A^.B"). 

35 9. A Base Station receiver according to one of the preceding claims, characterized in that said reduction-complexity 
control means (17, 18, 19) being arranged for minimizing for all the K users ti>e absolute value of the difference 
between the parameters of said temporal and spatial whitened channel model H^(16) and corresponding reduced 
set of re-estimated parameters (A^, B H ), in so determining the number r k of said cascades (Sj, f 1f s n t r ) for 
each of the K users. k 

40 

10. A Base Station receiver according to any claims from 6 to 9, characterized in that it further includes means (20) 
for promoting spatial and temporal re-correlation, also termed inverse-whitening means, acting on said reduced 
set of re-estimated parameters (A^, B^) in order to reconstruct said channel model estimated matrix in a form 
based on said cascaded spatial and temporal FIR filters (A^ B for each of the K users. 

45 Hk 

11. A Base Station receiver according to claim 10, characterized in that said means (20) for promoting spatial and 
temporal re-correlation being arranged for multiplicating said reduced set of re-estimated parameters (A . B by 
a second spatial factor R H/2 which is the inverse of the first one R^ 2 - and by a second temporal factor (X^J-^ 2 
which is the inverse of thetirstone (Xpc")** 2 , obtaining reduced se?of re-estimated and inverse-whitened param- 

50 eters A k , B* for each of the K users, where A^ is a matrix including said r k beamformers and B^ is a matrix including 

said r k temporal FIR. 

12. A Base Station receiver according to claim 11, characterized in that said multi-channel noise and interference 
estimating means (22) being arranged for re-estimating the space covariance matrix (R^j of the sample collection 

55 of additive ambient noise and inter-cell interference among the various antennas, by subtracting from the base 

band digital samples (Y) of the signals received on the array the residual of the channel re-estimation with the 
reduced set of re-estimated and inverse-whitened parameters A^, B*, obtaining a re-estimated space covariance 
matrix ^ stored in the noise channel memory (21) of the receiver. 
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13. A Base Station receiver according to claim 12, characterized in that it further includes a spatial whithening filter 
(SWF) of the base band digital samples Y of the signals received on the array, arranged to multiplicate the signal 
Y by a third spatial factorR^ tf/7 ( obtaining a spatial-whitened. 

5 1 4. A Base Station receiver according to claim 1 3, characterized in that said means for promoting spatial uncorrelation 
(23) being further arranged for multiplicating said matrix A k by said third spatial factor fx"*, obtaining r k spatial re- 
whitened beamformersaVi .....a*., . 

15. A Base Station receiver according to one of the preceding claims, characterized in that said convolving means 
10 (11 b) including said re-estimated parameters of the channel model maintaining the same structure of r k cascades 

(s 1t ft,..., s n t r ) for each of the K users. 

16. A Base Station receiver according to one of the preceding claims, characterized in that it includes means for 
despreading said base band digital samples of data transmitted from all the K users in case the signal were trans- 

15 mitted multiplied with orthogonal spreading sequences associated to the users one to one, to perform the so-called 

CDMA technique, the despreaded symbols being said /(sequences (sF,....*^..., v*)) of the K users spatially and 
temporally matched to the channel. 

17. A Base Station receiver according to claim 16, characterized in that the re-estimated parameters (s A1 , s kJ ) 
20 of said temporal FIR equivalent to as many signatures are obtained by convolving the stored temporal re-estimated 

parameters (b^, ...,b fcf ) by said spreading sequence (c k ). 

18. A Base Station receiver according to claim 17, characterized in that said multi-user detecting means (10) being 
arranged for calculating a correlating matrix (R) constructed by Kr x Kr spatial sub-matrix (R^ containing the 

25 correlations between the Kr spatial beamformers and Kr x Krtemporal sub-matrix (R^) containing the cor- 

relations between the temporal re-estimated parameters (s k } , s* r ), and for inverting said correlating matrix (R). 
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